JLAB-THY-99-06 
UND-HEP-99-BIG 03 
hep-ph/9903258 



Pauli Interference in the 't Hooft Model: 
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Abstract 



^ ■ Pauli Interference in decays of heavy flavor mesons, a genuinely nonleptonic non- 

perturbative effect, is considered in the 't Hooft model. Analytic summation of the 
exclusive decay widths yields the same expression as obtained in the OPE-based ap- 
proach. Threshold-related violations of local duality in PI are found to be suppressed. 
A novel case is identified where the OPE effect is completely saturated by a single 
final state. 



The dynamics driving heavy flavor transitions are concisely expressed on the quark 
level in terms of CKM parameters and masses. Yet their impact on the observable 
weak decays of heavy flavor hadrons Hq is obscured by nonperturbative features of 
strong interactions. The principal elements of the model-independent treatment of 
inclusive decay probabilities based on Wilson's operator product expansion (OPE) 
U, were laid down almost 15 years ago @, ||]. Significant progress has later been 
achieved in developing this technique and establishing theoretical control over non- 
leptonic, semileptonic and radiative channels of practical interest (for a recent review, 
see Refs. [|J). The nonperturbative corrections which would scale like 1/mg are ab- 
sent from the truly inclusive widths; the leading nonperturbative corrections arise 
at order l/rnq and are "flavor-independent" corrections insensitive to the flavor of 
the spectator (s) [|5], EJ. Effects sensitive to the flavor of the spectator (s) emerge at 
the 1/rriQ level [|[ [7|, |5|. They are conventionally called Weak Annihilation (WA) 
in mesons ||, Weak Scattering (WS) in baryons and Pauli Interference (PI) |TD[ 
in both systems. OPE relates their contributions to the expectation values of local 
four-quark operators. 

The OPE is typically constructed as an asymptotic power expansion in Euclidean 
space, and neglects the terms ~ exp (— v^FMqcd) suppressed for high momentum 
scales \fQ* ' . Yet in Minkowski space these can lead to terms which oscillate and have 
only power suppression. For heavy quark decays the practical OPE a priori does not 
forbid, for instance, dependence of the type sin (juqAqcd) / ^Iq 0, 0- ^° wm ch 
degree they are suppressed by powers of 1/rriQ depends on details of the strong forces 
and the specifics of the process under study. This problem is behind the violation 
of local quark-hadron duality. Since WA and PI represent power suppressed effects, 
one might expect numerically larger corrections to duality at finite values of niq than 
for the leading terms in the OPE. Even a more radical concern had been expressed 
in the past that the OPE is not applicable to nonleptonic heavy flavor decays at all, 
even at arbitrary large heavy quark mass. 

Implementation of duality in nonleptonic decays is not always obvious at first 
glance, in particular for PI. The underlying quark diagram for PI in B~ mesons is 
shown in Fig. 1. The u antiquark produced in the decay must be slow to interfere with 
the valence u. Tracing the color flow one is then faced with a dichotomy. The large 
momentum flows through the diquark loop (cd) which represents the "hard core" of 
the process. In practical OPE one effectively replaces the propagation of this diquark 
by a nearly free di-fermion loop; the absorptive part is evaluated - at least to leading 
order - via the production of free quarks. Such identification cannot be justified by 
itself even at arbitrarily large momentum transfer, since the spectrum of QCD does 
not contain colored states. 

Alternatively, one can combine a hard quark q from the loop with the slow valence 
antiquark to arrive at a color singlet configuration possibly dual to the hadronic 
spectrum. That comes at a price, though: this qq sp system is no longer a hard one. 
For at least in a partonic picture with p sp ~ m sp — > the qq sp invariant mass vanishes 
irrespective of m b . 
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Figure 1: Quark diagram for PI in B decays. The weak vertices are broken to show 
the color flow for the leading- N c contribution. 



An intriguing example that interference might be delicate was noted in Ref. [[13 
where an inclusive and exclusive treatment of B decays seemed to yield different 
results for the interference of two separate color operators. 

A dedicated consideration suggests that such apparent puzzles can be resolved 
without invoking new paradigms or assuming a higher onset of duality, growing with 
iV c . Here we will present an explicit analysis of PI in a soluble field theory, the 't Hooft 
model |TJJ] which is QCD in 1+1 dimensions in the limit N c —>-oo. We have found that 



the OPE expression matches the result obtained when summing over the hadronic 
resonances; the onset of duality is largely independent of the details. Moreover, 
the parton-deduced OPE expression appears to be exact in the chiral limit when all 
involved quarks except Q are massless. In this case the considered effect is saturated 
by a single final state. From this perspective it is complementary to the classical 
small velocity limit in semileptonic decays || where a similar saturation holds in the 
actual D = 4 QCD for heavy enough b and c quarks, up to power corrections. 

The 't Hooft model has already been employed recently for studying inclusive 
heavy flavor decays, both analytically [|T5[] and numerically The large- iV c limit 



harnessed in the analyses, however, washed out the difference between semileptonic 
and nonleptonic decays. In the case of weak annihilation, additionally, it equates the 
validity of the OPE for the heavy flavor width with its applicability to the vacuum 
current-current correlators, and thus can be thought not to be representative for the 
problems encountered in the nonleptonic decays. PI is a genuinely nonleptonic effect. 
Invoking the spectator quark makes this a leading- N c effect and allows a consistent 
evaluation in the framework of the 't Hooft eigenstate problem. Presentation of full 
analysis and details of the heavy quark expansion in the model is left for a separate 
dedicated publication. 

The 't Hooft model and heavy quark decays 

The £ t Hooft model - QCD in 1+1 dimensions with iV c — >oo - has been described in 
many papers [jrj|, [IT]]. While the Lagrangian has the usual form 
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the coupling g s has dimension of mass making the theory superrenormalizable. 
has dimension of mass as in D = 4 in the adopted normalization; however the fermion 
fields i/i(x) carry dimension of m 1 / 2 . In the large N c limit strong interaction effects 
are driven by the parameter 

playing the role of the nonperturbative scale Aqcd- Additional benefit of the D = 2 
theory is that the observables can be expressed in terms of bare masses and couplings, 
and that the perturbative corrections for hard quantities generate power series. More 
relevant details and discussions can be found in Refs. [TT5]. [TJ3 . 



In the limit N c ^oo the spectrum of 1+1 QCD consists of mesonic quark- ant iquark 
bound states which are stable under strong interactions. The meson masses are given 
by eigenvalues of the 't Hooft equation 



M n V„(x) 



m\ — (3 2 ml — f3 2 



x 1 — x 



y n {x) - f 1 dy , (3) 

Jo [y — xy 



where m 12 are the masses of the quark constituents, and the integral is understood 
in the principal value prescription. The solutions are the light-cone wave functions 
<p n (x). 

Polarization operator of the vector dj^u current 

IV(g 2 ) = -U(q 2 ) (q 2 5^-q,q u ) , p{q 2 ) = --\mll{q 2 ) (4) 
(assuming that m u = irid and the current is strictly conserved) has a resonance form 

n(? 2 ) = - £ -J^— , p(q 2 ) = vr £ f 2 5 (q 2 - Ml) , (5) 

n y n n 

where the decay constant /„ of a particular meson n is given by 

fn = \ — [ dxip n (x) . (6) 



7T JO 

In D = 2 the exact correlator of vector currents for massless quarks is very simple: 

I% 2 ) = ^, p(q 2 ) = N c 5(q 2 ) . (7) 



In the 't Hooft model at m u d = one has fo = y N c /tt and M = 0; for all excitations 
/„ = 0. With nonzero quark masses the spectral density shifts upward, to the mass 
scale ~/3m tt (i or m 2 d , and a high-energy tail in p appears ~N c (m 2 + m 2 d )/q A |TJ 
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We choose the weak decay interaction of the current-current form. Since in 1+1 
dimensions the axial current is related to the vector one, = e^ u J^, we simply 
consider the VxV interaction, and include both charge- and neutral-current vertices: 

G - - 

£ WC ak = —j= {a^crff.b) (d'fu) + a 2 {d lfl b) {cYu)) + H.c. . (8) 

While G is an analogue of the Fermi constant, it is dimensionless here. In the example 
of the parton width the interference term oc2aiCi2 is non-leading in 1/N C : 

r quark = -— N c (a^ uark + ^rf ark + 2a!a 2 r^ ark ) , (9) 

where T^ ark are ~0(N®), while Y 1 ^^ ^^(l/Nc). Yet also the interference term can 
be leading in 1/N C |19[ - if the contributions are non-leading in l/m^. Specifically, 
considering B~ decays affected by PI, see Figs. 1 and 2, we find Y\l ~ 0{N®)\ the 
three terms Y\, Y 2 and Y\ 2 l are all of the same order in 1/N C . This allows one to retain 
only the leading- N c decay amplitudes. 




a b 



Figure 2: Large-iV c decay amplitudes induced by charge-current (a) and neu- 
tral-current (b) terms in the weak decay Lagrangian. 

There is something else we can read off from Fig. 2. The two quark-antiquark 
clusters in the final state are produced with quite different characteristics: while the 
upper one is produced by a "pointlike" source, the lower is not, since it combines a 
decay quark with the valence quark; we will refer to the latter as "multiperipheral" 
production. On general grounds one expects the typical mass distribution for hadrons 
from the upper vertex to be considerably different from that from the lower one. 
Normally, the former is much wider spanning over the whole accessible range up 
to m&. However, for the special case of vector-like interactions in D = 2 the bulk 
of pointlike-produced hadrons has smaller masses comparable to the corresponding 
quark masses. Yet interference between the two diagrams in Fig. 2 requires the du 
and cu clusters to have the same mass in both cases. This explains why PI is power 
suppressed - a property that emerges automatically in the OPE. 
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We turn now to the 't Hooft model to transform these qualitative observations 
into specific expressions. Fig. 1 generates - to leading order in the OPE - the operator 



G 2 

-aia 2 — 
A 



m 2 b / 



2m c m d 



mi 



(bu) (ub) + {bi^^u) {ui~j b b) 



K 



1 



lm c + m d ) 



1 



[m r 



m d y 



mi 



mi 



1/2 



(10) 



This contribution remains finite for m c = m d = m u = and can conveniently be 
analyzed in this chiral limit. The B~ expectation value of the operators in Eq. (|i~0|) 
has a iV c -favorable color structure and is, therefore given by vacuum factorization: 



2M B 

thus yielding in the chiral limit 



(B-\(b ltll5 u) (« 7AtT5 6)|S-> 



•.pi 



q2 

-2a x a 2 — f 2 B M B . 



(12) 



We note that T PI approaches a constant for m b ^oo since fs oc 1/ y/rrib. It constitutes 
a 1/mq nonperturbative correction to the decay width, in contrast to 1/ttiq in four 
dimension. This follows from different canonical dimension of quark fields |T5| . 

There are no 1/ m;, corrections to this result. This is a peculiarity of two dimensions 
where the absorptive part of the (di) quark loop in Fig. 1 scales as the momentum to 
the zeroth power and, thus does not depend on whether quark or meson momentum 
flows through it. The corrections to the Wilson coefficient or other higher-order 
operators can induce only terms suppressed by at least two powers of l/m b , if m c>u ^ 
m b . 

Let us now consider the decays in terms of hadrons. To leading order in 1/N C 
the final states are pairs of mesons generically referred to as D® and 7r~. The partial 
decay width B 



Dln n takes the form 



G 2 



kn 



8M 2 \p\ 



a\\A k B n \ 2 + a\\A n B k \ 2 + 2a x a 2 Re A k A* n B* k B n 



(13) 



where A and B schematically denote the "multiperipheral" transition amplitudes and 
the "pointlike" meson creation amplitudes, respectively: 



Ai 



B„ 



(n\J,\0) 



(14) 



With p denoting the rest-frame momentum of the final state mesons, the PI width is 
then given by the sum 



G 



2a i«2 ^rjf E i^T A k A* n B* k B n 



8M 2 B W\ 



(15) 
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Eqs. ( |T3| , |i"5| ) reflect the peculiar two-body phase space ocl/|p | in two dimensions. 

The quantitative matching between the OPE-based calculation and the hadronic 
saturation of the interference width is most easily seen when all final state quarks 
u, d, c are massless. While not affecting the OPE analysis, this limit significantly 
simplifies the expressions for the individual hadronic amplitudes, as explained in 
Refs. fl5| , |18|j . In the case at hand only n = survives for the decay amplitude ~ a x 
(Fig. 2a) and k = for the amplitude ~ a 2 (Fig. 2b). The interference then resides 
in the single final state containing the lowest lying massless D° and 7r~. Moreover, 
the transition amplitudes between the two mesons take a particularly simple form at 



q = [15| in terms of their 't Hooft wavef unctions: 



1 r 1 

%^{ k \^uJu\H Q ) = -q z dx <p k (x)<p B (x) , (16) 

with the 'pointlike' amplitude given by £> ~ ze^/o-Pj yielding 



,pi 9 G 2 M B N C 
= —2aia 2 



47T 







dx ip B (x) 



2 = -2a x a 2 ^ f B M B , (17) 



where we have used the fact that <po(x) = 1 for massless quarks. The minus sign 
emerges since the direction of the vector playing the role of q is opposite for the two 
interfering amplitudes. 

Thus, the OPE prediction Eq. (||) is exactly reproduced. Apparently, there is no 
violation of local duality at all for PI when m u = md = m c = ! It could be anticipated, 
for in this limit the only threshold in T PI occurs at mb = 0. Then the OPE series can 
have the same convergence properties in Minkowskian as in Euclidean space |TT . 



With m Uj( i, c 7^ the situation becomes more complex since more production thresh- 
olds arise. Those exhibit singularities due to the singular two-body phase space in 
D = 2. As explained in detail elsewhere, the predictions of the practical OPE must 
be compared to the decay probabilities where the threshold singularities are aver- 
aged, and for heavy quark decays this procedure is naturally done by smearing over 
the interval of the heavy quark mass, exceeding the distance between the successive 
principal thresholds O, |I3fl . Since this phase space factor is integrable, the threshold 
spikes do not affect the width smeared over a mass interval Ami, ~ 1/^6- 

The calculation uses the detailed kinematic duality between the partonic and 
hadronic probabilities. Details of the derivation for the massive case can be found in 
Ref. [pT8|| . Here we only sketch the basic steps. First, one takes into account that for 
the bulk of the decays the masses of both 'pointlike' i and 'multiperipheral' j mesons 
are small compared to m^. 

M 2 £ m 2 q , f3m q ; M] ^ f3m b , (18) 

where m q generically denotes the final state quark masses. These features anticipated 
from the quark-level consideration has been demonstrated directly in the 't Hooft 
model in Refs. |15|, where their accuracy was also quantified. It is sufficient for 
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calculating T up to terms 1/m 2 . Then we can expand momentum \p\ and transition 
amplitudes around the massless kinematics: 



\P\ 



M 



1 + 2 



B 



M 2 k + M 2 
Mn 



+ 



M 2 

A(o) + -£ 



mi 



m\ 



dAt 



dq 2 



(19) 



The set of sum rules derived in Ref. |To] allows to calculate the sum of the decay 
probabilities and the leading corrections associated with nonzero meson masses; the 
latter cancel the increase in l/\p\ in Eq. (|15|). 

For example, neglecting deviation of q 2 from in calculating the 'multiperipheral' 
amplitudes, we get for the two interfering amplitudes [O, ITS 




{M 2 B - M 2 k ) 



dy Vn{y) 



(M 2 B -M 2 ) dx<p B {x)ip n {x)- dy<p k {y). (20) 
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Extending summation over all k and n and using the completeness of the 't Hooft 
eigenstates, we obtain 



-ipi 



pi 

kn 



-2aid2 



G 2 N C 
4nM B 



[ dxif B (x) 
Jo 



9 02 f 

-2a x a 2 — J 



(21) 



This expression is valid up to 0(l/mf) corrections. The leftover effect of the slope of 
the transition formfactor is quadratic in m q /mb [ [Tj| . 

Thus, T PI agrees with the expression given by the free quark loop to the accuracy 
suggested by the OPE. 

For small, yet nonvanishing final state quark masses one can estimate the limita- 
tion in local duality for PI as it relates to thresholds. Close to the opening of such 
a threshold the hadronic width is not exactly reproduced by the OPE. These effects 
are characterized by the nature of the threshold singularity, the scaling of its residue 
with mj and the distance between principal thresholds. 

It turns out that these effects essentially depend on the relation between the final 
state masses. The strongest violation of local duality occurs when one of the mesons 
is a low excitation while the other has a large mass close to m&. 

The case with m c y^0, m u = md = is somewhat special. For one of the interfering 
decay amplitudes vanishes right at threshold and T PI undergoes a finite jump there 
(details can be found in Ref. fT8|j): 



5T^q\ oc —2a 1 a 2 fBm l 



G 



02 



7T 2 /3 2 



M 



9/2 



9{M B -M k 



M k+l -M k 



2M 



(22) 



B 



A more elaborate analysis suggests that the relative sign of the two interfering am- 
plitudes alternates for successive thresholds, and we arrive at the following ansatz: 



ST 



PL 

osc 



oc — 2a x a 2 fern, 



c mT 



J2(-l) k 0(M B -7f(3Vk) . 



(23) 
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The oscillating term in PI thus gets damped by at least l/m\. 

For m U ( i nonzero, the picture changes essentially in two respects: neither decay 
amplitude vanishes at threshold, however the phase space factor softens to [2m n (M B — 
M t hr)]~ 1//2 for the decays B — > D^> + 7r at m n ^ 0. One then finds 



AS oc % 2a ia2 EH)* "^"^ , M£ * . (24) 



2 ■ m b r v ft^E 



With all final-state masses non-zero, many additional thresholds open up and the 
dynamical landscape becomes in general rather complex. Yet one observes that the 
contributions from the final states where both meson masses constitute a finite frac- 
tion of rrtb are suppressed by even higher powers of 1/ m& 



Conclusions 

Let us first describe the hadronization picture suggested by the large-m^ analysis. 
We observe that duality applies to the quark-antiquark pairs where each energetic 
quark (c or d) is combined with the wee spectator antiquark or slow u produced in 
the weak vertex. The completeness of the hadronic states - or, in other words, the 
duality between the parton-level and mesonic states - is achieved already for a single 
fast moving decay quark when it picks up a slow spectator. The 'hardness' of these 
processes determining the applicability of the quasifree approximation, depends on 
the energy of the fast quark rather than on the invariant mass of the pair. For actual 
hadrons, the slow spectator (s) compensate the color of the initial static heavy quark. 
This initial distribution of the color field singles out the rest frame and introduces 
the new scale, energy vs. invariant mass. 

In the full graph for the meson decay which would explicitly include propagation 
of the spectator, there is always a color mate for any quark in the "partonic" part 
of the diagram. Therefore, considering the colored diquark loop as nearly free yields 
the correct result for the properly formulated problem. 

In this note we have analyzed PI as a spectator-dependent nonleptonic effect for 
heavy flavor decays within the 't Hooft model. We have found the OPE expression 
to match the result obtained when summing over exclusive final states at least up 
to terms (AqcdA^&) 3 T_b which is the accuracy we adopted. This is the first anal- 
ysis where the peculiarities of nonleptonic decays widths of heavy flavors emerge 
compared to their semileptonic decays or processes like e + e~ annihilation. We did 
not encounter specific problems for applying the OPE, in accord with the general 
arguments justifying it for all types of fully inclusive decay widths [0, 0, [T5|, [L8|]. 

We estimated the effects of violation of local duality in PI in the heavy quark 
limit. It turned out to be strongly suppressed scaling like l/m\. A novel case was 
identified where local parton-hadron duality is exact. Although applicable only in two 
dimension for large N c , it complements the classical case of small velocity semileptonic 
b — > c decays 0| which was instrumental in developing the modern applications of 
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the heavy quark expansion. 
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